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Abstract. Photometry of the Q2237+0305gravitational lens in VRI spectral bands with the 1.5-m telescope of 
the high-altitude Maidanak observatory in 1995-2000 is pre- sented. Monitoring of Q2237+0305 in July-October 
2000, made at nearly daily basis, did not reveal rapid (night-to-night and intranight) variations of brightness of 
the components during this time period. Rather slow changes of magnitudes of the components were observed, 
such as 0.08m fading of B and C components and 0.05m brightening of D in R band during July 23 - October 
7, 2000. By good luck three nights of observation in 1999 were almost at the time of the strong brightness peak 
of image C, and approximately in the middle of the ascending slope of the image A brightness peak. The C 
component was the most blue one in the system in 1998 and 1999, having changed its (V-I) color from 0.56m to 
0.12m since August 1997, while its brightness increased almost 1.2m during this time period. The A component 
behaved similarly between August 1998 and August 2000, having become 0.47m brighter in R, and at the same 
time, 0.15m bluer. A correlation between the color variations and variations of magnitudes of the components 
is demonstrated to be significant and reaches 0.75, with a regression line slope of 0.33. A color (V-I) vrs color 
(V-R) plot shows the components settled in a cluster, stretched along a line with a slope of 1.31. Both slopes 
are noticeably smaller than those expected if a standard galactic interstellar reddening law were responsible for 
the differences between the colors of images and their variations over time. We attribute the brightness and color 
changes to microlensing of the quasar's structure, which we conclude is more compact at shorter wavelengths, as 
predicted by most quasar models featuring an energizing central source. 

Key words, cosmology; gravitational lensing - galaxies: quasars: individual: QSO 2237-1-0305 - methods: observa- 
tional - techniques: image processing 



1. Introduction 

The Q2237+0305 gravitational lens (the Einstein Cross) is 
one of the most impressive manifestations of the gravita- 
tional lensing phenomenon - four images of the same high- 
redshift quasar [z — 1.695) are arranged almost symmetri- 
cally around the lensing galaxy nucleus (z = 0.039) within 
a circle of approximately 2" diameter. The Q2237+0305 
system is an excellent target to study microlensing events, 
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because the light beams, corresponding to the 4 lensed 
quasar's images, pass through the interior, heavily pop- 
ulated part of the lensing galaxy and thus, have a high 
probability to intersect a significant mass of microlens- 
ing stars as they pass through the inner disc (Kayser & 
Refsdal 1989). 

The system has been intensively examined since 1987, 
when the first published measurements of magnitudes of 
the individual lensed quasar components in g,r and i 
Gunn filters were made by Yee H1988|l . The first attempt to 
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build the light curves of the four quasar components was 
made in 1991 by Corrigan et al. (,199 L) ■ They brought to- 
gether all the available Q2237 images of sufficient quality, 
taken with different telescopes and in a variety of pass- 
bands, - Mould B, V and R and Gunn g, r and i, - and 
reprocessed with a single algorithm. Having used the mul- 
ticolor photometry data for 33 normal stars, whose (B—V) 
colors ranged from -0.3 to 1.5, they calculated the relevant 
color equations, which allowed them to reduce all the ob- 
servations to a single passband. Their r Gunn and B light 
curves cover the time period from September 1986 up to 
December 1989, and include the first microlensing event 
observed in August 1988 by Irwin et al. H1989|l . 

A further attempt to use all the available observational 
data for Q2237-I-0305 was made in 1994 by construct- 
ing "differential" light curves, which were argued to be 
free from the effects of different spectral bands, technique 
of zero-pointing, and the quasar intrinsic brightness vari- 
ations, (Houde et al. 11994(1 . In addition to the data of 
1986-89, contained in (Corrigan et al. 1991J, other results 
were used, taken in 1990 and 1991 by Crane et al. (^91j, 
Racine iflM^ . Rix et al. iflM^ . and Houde et al. (T994). 

The first program of regular photometric monitoring 
of Q2237-h0305 was started in 1990 at the Nordic Optical 
Telescope, (0stensen et al. ll996f) . A large number of mea- 
surements of the four quasar components in V, R and / 
spectral bands during five years were obtained, which per- 
mitted construction, with the use of Corrigan et al. H1991|l 
zero-pointing, of the historic light curves, covering 9 years 
of observations. 

No regular multi-filter monitoring of Q2237-I-G3G5 is 
reported between 1996 and 1999, excepting our results 
of VRI photometry for three nights on 17-19 September 
1995, (Vakulik et al. I1997|l . and the similar results by 
Burud et al. I|1998|l . obtained with the Nordic Optical 
Telescope for a close epoch, 10-11 October 1995. No night- 
to-night or intranight brightness variations of the four 
components have been found for these time periods, while 
a noticeable change in the component B color as compared 
to the observations by Yee in 1987 H1988|l has been re- 
ported in both works. A short time-scale monitoring with 
the CFHT in June 14-16, 1992 should be mentioned here, 
(Gumming & De R.obertis. I1995|l . which also did not re- 
veal any photometric variations in R and / bands during 
a three-day period. 

In 1999 and 2000, the results of VRI photometry in 
1997 and 1998 with the Maidanak 1.5-m telescope were 
pubHshed, (Bhokh et al. 1999 ) and (Dudinov et al.l^OOOa). 
Recently, the superb results of a detailed long-term moni- 
toring, obtained within the OGLE program from 4 August 
1997 to 5 November 2000 has become publicly accessi- 
ble, and partly presented in two papers of Wozniak et al. 
iPinOa,) and ( 2000b). And in 2002, the results of moni- 
toring of Q2237-I-0305 by GLITP collaboration appeared, 
which cover the 4-month period October 1999 - February 
2000 (Alcalde et al.|2002). The most recent publication of 
the results of low-resolution observations with the 3.5-m 
telescope at the Apache Point Observatory should be also 



mentioned here, (Schmidt et al. I2()()2|l . which have given 
the Gunn r lightcurves of A and B components for 73 
dates between July 1995 and January 1998. In spite of a 
rather low photometric accuracy, - an error bar of 0.1™ to 
0.2™ is reported, - the data are of value first of all because 
they include the brightness peak of A component in 1996. 
As far as we know, no more data about this event have 
been ever published, though reported in private commu- 
nication, (e.g. R. Ostensen). 

The brightness records taken in broad-band filters are 
a good starting point for theorists to estimate the size 
of the quasar radiating region in the visual (ultraviolet 
rest frame) continuum, and to determine the range for 
microlens masses responsible for the observed brightness 
variations. Both the statistical analysis of long-term mon- 
itoring data, and simulation of the isolated microlens- 
ing peaks have been applied to calculate these values, 
(Nadeau et al. lBMl Lewis & Irwin HlMl Refsdal & Stabell 
Eni Webster et al. ITMTI Wyithe et al. ''iOOOa', Wyithe 
et al. I2000bl Wyithe et al.'2002, Yonehara 2000), having 
given the estimates of the quasar dimension, - lO^^cm to 
lO-'^^CTO in the optical continuum, and a great variety of mi- 
crolens masses ranging from O.OOO6M0 < M < O.OO6M0, 
(Nadeau et al. [TMT|l . to O.IMq, f Wambseanss [TMT|) . all 
indicating however, that microlensing events in the system 
are mostly caused by subsolar-mass objects. 

These available models confront a problem which prob- 
ably indicates a breakdown of the simple assumption that 
the luminous quasar is a uniformly bright accretion disc. 
The dilemma is that the short time scale of the observed 
events like the C image peak in July 1999 and the A bright- 
ness peak in November 1999 occur on such short time 
scales that small accretion disc diameters are implied; such 
small bright accretion discs would have occasional strong 
brightness peaks of several magnitudes that are never ob- 
served. This probably tells us that a more complex model 
where the time scale of the brightness peaks is related to 
a quasar structure crossing time, not the crossing time of 
the entire quasar luminosity. We expect to apply, in a sub- 
sequent report, the Schild and Vakulik (12003') double ring 
quasar model that successfully models the long history of 
Q0957+561 microlensing observations. 

No systematic multicolor photometric measurements 
existed until the monitoring program with the NOT was 
started in 1990 (Ostensen et al. I1996|l . Meanwhile, a sus- 
picion was expressed by Corrigan et al. in 1991, and in 
1992 by Rix et al. independently (Corrigan et al. I1991L 
Rix et al. ll992|l . that the color indices of the components 
might have changed since the first three-color observations 
by Yee, (1988^. R was a very important statement, since 
reducing different datasets to a single light curve, (e.g. 
Houde et al.,'1994!), as well as determining the extinction 
law in the Q2237+0305 lensing galaxy, (Yee JJi_88, Nadeau 
et al. 1991, Falco et al. I1999|l . are substantially based on 
the assumptions, that "all four components have identical 
intrinsic color indices", and "their observed color differ- 
ences are due to different degrees of interstellar extinc- 
tion and reddening by the same extinction law" , (Houde 
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et al. I1994|l . and "the magnification is wavelength inde- 
pendent... and time independent", (Falco et al. 1999 ) . In 
particular, Falco et al. IflD^gt measured the value of Ry 
for Q2237+0305 to be equal 5.3 and came to a conclusion 
about great differences in the extinction laws for lensing 
galaxies from a sample consisting of 23 gravitational lens 
systems. However, as it can be seen from Fig. 4 in Falco 
et al. (199.9J, the differences may be significant at wave- 
lengths shorter than 550nm, while at larger wavelengths 
the difference between the extinction curves does not ex- 
ceed the error bars. 

In discussing the results of the five-year VRI moni- 
toring of Q2237-I-0305, 0stensen et al. H1996|l did not an- 
alyze, however, any color changes in the system, having 
noted only " very nearly equal" colors for the components 
A and B, as well as roughly equal colors of C and D, with 
the extinction difference between the pairs of 0.6™ in V 
band, provided the extinction law follows A^^, according 
to Houde et al. fTOT). Meanwhile, Vakuhk et al. l(TM7|l 
and Burud et al. (,1998j) reported that the B component 
became the most blue one in the system in 1995, as com- 
pared to observations by Yee in 1987, fYee ll98H|l . 

The next step in determining colors and color changes 
in the Q2237+0305 system was made in VRI observa- 
tions with the Maidanak 1.5-m telescope in 1997-1998, 
presented in Dudinov et al. .2000a. and Dudinov et al. 
I2000bl Variations of colors were argued to be significant, 
and a tendency of the components to become bluer as their 
brightness increased has been demonstrated with the use 
of all available multicolor data. Unfortunately, the remark- 
able monitoring by Wozniak et al. I2000al and Wozniak et 
al. l2000bl has been made only in V band, and thus can not 
be used to investigate the color changes, while the most 
recent data of the GLITP collaboration have been taken in 
V and R filters for a campaign of 4 months only (Alcalde 
et aLEOnS). 

By this time, a great amount of observations of 
Q2237-I-0305 in spectral ranges other than visual contin- 
uum exists, - VLA observations at 20cm and 3.6cm (Falco 
et al. 1996), observations in the near and mid-IR (Nadeau 
et al. 119911 Agol et al. I2U0U|I . and in the quasar emission 
lines (Fitte & Adam lTUMl Racine lTM^ De Robertis & Yee 
Lewis et ai nMSl Saust lTUMji . The observed mag- 
nitudes of the components have been found to be almost 
unaffected by microlensings in these spectral ranges, which 
indicates that much larger quasar features radiate in IR 
and in the radio, as well as in the emission lines, as com- 
pared to the optical continuum. The recent detection of an 
arc of C III] emission, connecting A, D and B components 
(Mediavilla et al. I1998fl . should be regarded as a visual 
proof of the extended emission line region of the source. 
Because of the low sensitivity of a large source bright- 
ness to microlensing, the brightness ratios for the compo- 
nents, measured in these spectral ranges, were used to test 
the validity of a great variety of the existing macrolensing 
models, listed by Wyithe et al. ^UU^. 

Observations in UV with the HST (Blanton et al. 1998 ) 
and the recent X-ray imaging of Q2237-I-0305 with the 



Chandra X-ray Observatory (Dai et al. I2003 >) should be 
also mentioned here, which provided, in particular, highly 
accurate relative coordinates of the components (Blanton 
et al. T998) and the upper limits for the physical size and 
brightness of the Broad-Line Region producing Ly-a emis- 
sion, (Dai et al. I2()()3(l . Also, the Chandra data permitted 
calculation of the time delay between the A and B com- 
ponents of 2.7 hours. 

2. Observations 

Our observations were carried out with the 1.5-m AZT- 
22 telescope of the high-altitude Maidanak observatory, 
(Central Asia, Republic of Uzbekistan), known for its su- 
perb seeing conditions and a large number of cloudless 
nights, (Ehgamberdiev et al. I20()0|l . Because of technical 
reasons, we had to use three different CCD cameras in our 
observations, Pictor-416 camera in 1995, Pictor-416 and 
TI 800 X 800 cameras in 1997 and 1998, and ST-7 camera 
in 1999 and 2000. And because of technical reasons again, 
both f/8 and f/16 focal lengths were used in observations. 
The LN-cooled TI 800 x 800 camera, with pixel size of 
15/i, kindly provided by Prof.D.Turnshek, unfortunately 
revealed some peculiarities, caused by the charge trans- 
fer inefficiency, that is characteristic for the CCD's of this 
generation, (Turnshek et al. 11997(1 . In particular, notice- 
able stretching of stellar images in the direction of charge 
transfer is observed, as well as a dependence of the PSF 
upon coordinates at the chip plane. In addition, sensitivity 
irregularities of the chip can not be corrected satisfactorily, 
with the output of the flatfielding procedure dependent on 
the signal level. All these peculiarities reduced the actual 
accuracy of photometry, that is seen in Tables Eland El 

Unfortunately, a poor telescope tracking system 
spoiled the intrinsically good seeing of Maidanak site 
sometimes and did not permit use of exposures longer 
than 3 minutes. To provide sufficiently high accuracy of 
our photometry with such short exposures, we took images 
in series, consisting of 10 to 20 frames each. The frames 
were averaged before being subjected to photometric pro- 
cessing, while a comparison of photometry of individual 
frames enabled us to obtain an adequate estimate of the 
random error inherent in a particular series. 

Most of images has been taken in R band, - 31 dates 
in 1995-1999, (Tabled, plus 46 dates in 2000, (Table El 
- which were obtained almost at a daily basis during 2.5 
months. There is also photometry in V and / bands for 17 
dates in 1995-2000, (TableEJ- Some results have been pre- 
sented in our previous publications (Vakulik et al. 119971 
Bliokh et al. ITMl Dudinov et al. gOOOa, Dudinov et al. 
I2000b|l . We present here the results of all our observa- 
tions, including those which have been never published. 
In particular, the observations of July-October 2000 are 
presented, which have been undertaken to search for short- 
period (night-to-night) variations of brightness. The ap- 
pearance of the Einstein Cross at six epochs between 
October 1995 and August 2001 can be seen in Fig. 1, 
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which clearly demonstrates high photometric variability 
in the system. 

In addition to magnitudes of the components, the see- 
ing conditions are also presented in Tables 01 and [S] - the 
values of FWHM for particular nights, the scales and the 
CCD camera used. 

3. Photometric Reductions 

The difficulties inherent in accurate photometry of 
ground-based images of Q2237-I-0305, have been noted by 
many authors, (Burud et al. ITOMl Corrigan et al. [TMTl 
Vakulik et al. HWl Yee irUggjl . They are due mainly to 
its extremely compact spatial structure, with the wings 
of the quasar images overlapping even under good seeing 
conditions. Additional difficulties are due to the presence 
of a rather bright foreground lensing galaxy, with its steep 
radial brightness distribution. These are the main reasons 
for poor agreement of the results of different monitoring 
programs, and even for a noticeable discrepancy in pho- 
tometric results for the same data reduced with different 
algorithms, (Burud et al. Alcalde et al. In 

photometry of the data of 1995-1999, we used the method 
described in Vakuhk et al. [TM71 and Bhokh et al. 1999 , 
which is in general features similar to the double iterative 
PSF subtraction method, proposed by Yee (4988) . who 
was the first to present spatially resolved photometry of 
the system. In short, the method consists of the following. 

The PSF estimate is obtained from a reference star 
image, and is further superimposed upon each image com- 
ponent and the galaxy nucleus alternately, and then sub- 
tracted in such a way, that no depressions would appear 
in the residual brightness distribution. Such a procedure 
is repeated iteratively until a stable convergence of esti- 
mates of brightness and coordinates of the components 
is achieved. Then, according to the resulting estimates, 
the quasar components are subtracted, and the residual 
galaxy brightness distribution is smoothed with a rather 
broad median filter. After the resulting galaxy bright- 
ness distribution is subtracted from the initial image, re- 
moval of quasar components is repeated, followed again by 
smoothing the galaxy brightness distribution with a suc- 
cessively decreasing window. The iterative process stops 
when the width of the median filter becomes of order of 
the PSF width. 

In processing our data for the 2000 observing season, 
another method was applied, which used the known rel- 
ative coordinates of the components and an analytical 
model of the brightness distribution in the galaxy, repre- 
sented as a sum of three two-dimensional Gaussian func- 
tions. Before describing the algorithm, consider the basic 
principles of photometry for compact groups of star-like 
objects, that have been implemented in the known algo- 
rithms of other authors. 

Even in the images of Q2237-I-0305 taken with the 
Hubble Space Telescope, the quasar components are star- 
like and thus, in the isoplanacy region, the entire picture 
(photometric model of the system) can be represented as 



a sum of the PSF's r{x — Xk, y — yk), and the galaxy hght 
distribution g{x,y), and in the case of a sampled CCD 
image, may be written as: 

4 

^^hr{i - Xk,j - yk) +9{hj), (1) 

k=l 

where i and j are pixel numbers in x and y axes, chosen 
in parallel to the CCD lines and columns, respectively. 
The unknown parameters, - the coordinates of the compo- 
nents in the detector reference frame, Xk,yk, their relative 
brightnesses Ik, and the galaxy light distribution g{i,j), - 
are usually estimated from a requirement to minimize the 
difference between a model and the observed brightness 
distribution in the detected image according to some cri- 
terion, - e.g. the minimum of the sum of square residuals 
criterion: 

= - = (2) 

i 3 

Here F{i,j) is the brightness distribution in the detected 
image, and the set of unknown parameters is denoted as p 
for short. The estimate of the PSF can be obtained from 
the images of reference stars near the object. 

As was noted above, noticeable difficulties in photom- 
etry of Q2237-I-0305 components are caused by the fore- 
ground lensing galaxy, with its light distribution g{x,y) 
being unknown. In minimizing Eq. |21 or another one sim- 
ilar to it, - the galaxy brightness distribution is usually 
represented either analytically, - (e.g. Burud et al. I1998L 
Alcalde et al. 12002(1 . or its digital form g{i,j) is estimated, 
- e.g. the MCS algorithm (Magain et al. ITMH Burud et 

ai. rnm . 

To solve the problem, iterative algorithms are often 
used, which in fact approximately realize minimization of 
Eq.[51 and also permit to obtain the estimate of g{x, y) ei- 
ther analytically fTeuber ll993[ Ostensen et al.'19 96|l or in 
a digital form, (Yee 196HI Vakulik et al.|l997,). The result- 
ing analytic or numerical model can be treated further 
in photometry of Q2237-I-0305 components as a known 
function. Such an approach noticeably simplifies the solu- 
tion procedure, and provides good intrinsic convergence, 
(Corrigan et al. ITMTl Alcalde et al. i2002i Burud et al. 
.1 998) , but unfortunately, does not ensure the absence of 
systematic errors in estimating the magnitudes of the com- 
ponents caused by an inadequate galaxy model. 

A new image subtraction method proposed by Alard 
& Lupton H1999|l and successfully applied by Wozniak 
et al. fSODUi I2000b|l and by Alcalde et al. ipHiT^ in 
Q2237-)-0305 photometry, is seemingly free from this weak 
point. However, a comparison of photometry results for 
Q2237-h0305 published by the OGLE group and those ob- 
tained with other methods, reveals some systematics in 
the components magnitudes, that is probably caused by a 
bias of brightness estimates in their reference image. 

The PSF is usually represented either numerically, or 
as an analytic function. In this work, the following ap- 
proach was used. We transformed all the detected images 
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to the same axisymmetrical Gaussian PSF with a pre- 
assigned parameter using the inverse hnear filtration 
procedure: 



Table 1. Parameters of the photometric model of the lous- 
ing galaxy for Q2237-I-0305 system. 



W 



Fo{uJi,UJ„) ■ R{uJl,Un) 
f{uJi,LOn) 



(3) 



Here i^o('-^/i ^^n) is the Fourier transform of the initial im- 
age, F(i,j) is the transformed (standardized) image, and 
W is the inverse Fourier transform operator. A complex- 
valued inverse filter w{uji,LUn) = is composed 
from the Fourier transform of the initial PSF r{i,j). A 
function R{uji,uJn) = exp[—a'^{ujf + forms the 
Fourier spectrum of the standardized image with the 
Gaussian PSF for the given parameter as - To construct the 
inverse filter, a reference star about 64" south-west from 
the quasar, denoted as a star in Gorrigan et al. (1991| 
was used. 

In doing so, we did not try to noticeably increase the 
resolution in the initial images, and used a transforma- 
tion that is a linear one, and, in contrast to non-linear 
filtration methods, retains photometric accuracy. To ex- 
clude dependence of the resulting PSF on the signal-to- 
noise ratio in the Fourier spectrum of a specific image, we 
also did not use any optimizing algorithms of image re- 
construction, such as e.g. the well-known Wiener filtering. 

Since the restoring filter is normalized to unity at the 
zero spatial frequency, such a transformation retains the 
integral brightness of an image, and thus the estimates of 
the components' brightnesses can be made in the units of 
the reference star brightness. 

With such standardized images created, the sum in 
Eq^can be represented as 

4 

sil, j) =J2lk exp{-[(z - Xkf + {] - ykf]l2ol}, (4) 
fc=i 

where Us is an effective width of the resulting PSF. 

The distribution of light over the galaxy was repre- 
sented by a sum of three two-dimensional Gaussian func- 
tions: 

3 

= X! -^mexp{-[(i-a;g)cos<^m-|-(j-i/g)sin(^„]^/27]^- 



- [-{i ~ Xg)sin93„j + [i - yg) cos (^,„]^/2e^}, 



(5) 



where Xg, yg are coordinates of the galaxy center, are 
normalizing coefficients, rjm and Sm are parameters deter- 
mining the characteristic widths of the Gaussian profiles 
along the major and minor axes respectively, with their 
meaning understood from Eq. 5, and finally, defines 
the major axes orientation. Therefore, the photometric 
model of the system f{i,j,p) in EqElcan be represented 
as a sum of two constituents, s{i,j) and g{i,j), which de- 
scribe the quasar components (Eq^J), and a photometric 
model of the light distribution in the lensing galaxy, (Eq. 



m 


/ 






P.A.(°) 


1 


0.875 ±0.021 


0.264 ± 0.045 


0.206 ± 0.038 


57 ±3 


2 


0.090 ±0.011 


1.260 ± 0.075 


0.790 ± 0.041 


78 ±2 


3 


0.035 ± 0.002 


5.440 ± 0.530 


2.840 ±0.110 


58 ±4 



Table 2. Relative angular positions of Q2237+0305 
A,B,G,D components and the galaxy center (G) from ob- 
servations of 2000. 



Component 


Aa(") 


Asn 


A 


0.000 


0.000 


B 


-0.674 ± 0.003 


1.679 ±0.004 


C 


0.624 ± 0.005 


1.206 ±0.004 


D 


-0.867 ± 0.008 


0.513 ±0.003 


G 


-0.085 ±0.014 


0.939 ± 0.006 



5). A set of 26 unknown parameters denoted as p, con- 
sists of four pairs of coordinates x^, and coordinates of 
the galaxy center Xg^yg, normalizing multipliers /fc, Im, 
the parameters 77™, Srm and orientations of axes of the 
three Gaussian components of the galaxy model. 

To calculate the parameters of the galaxy photometric 
model, as well as the coordinates of the quasar compo- 
nents, a set consisting of 14 best quality images was se- 
lected that was obtained on September 2, 2000 in R filter 
under the atmospheric seeing of O'.'S and better. The im- 
ages were averaged and reduced, through the inverse linear 
filtration procedure described above, to the Gaussian PSF 
with = 0'.'34, (FWHM of O'.'S). 

The least-squares algorithm was used to calculate the 
brightnesses and coordinates of the components and the 
parameters of the galaxy photometric model from the con- 
dition expressed by Eq.[51 It should be noted, that in such 
a way we obtain parameters of the galaxy model, that is 
the result of the convolution of an actual galaxy light dis- 
tribution with the Gaussian PSF with the given (Js = 0'.'34. 
Since Gaussian functions were adopted both for the PSF 
and for the constituents of the galaxy model, the decon- 
volved galaxy model parameters can be easily calculated. 
Such deconvolved parameters are presented in Table ^ 

In TableEl the relative positions of the B, G, D compo- 
nents and the galaxy center in the equatorial coordinate 
system, calculated from the 14 selected images with the 
procedure described above are presented. Our coordinates 
agree within 0'.'015 with those obtained from the HST im- 
ages (Grane et al. lBMl Blanton et al. niM|l . 

In the subsequent photometric processing of all the 
available data, every image was reduced to a "standard" 
PSF, and the corresponding quasar image brightnesses 
were estimated by minimizing the function ([SJ, with the 
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parameters of the galaxy model and the relative coordi- 
nates of the components being fixed, according to Tables 
n and 121 The a star from Corrigan et al. H1991|l was used 
as a secondary photometric standard, with its magnitudes 
taken from this work. 

Photometry of the image sets taken during a single 
night does not show brightness variations that might be 
regarded as significant as compared to the photometry 
uncertainties. Therefore, the brightness estimates taken 
within a night were averaged, and the formally calculated 
error in the mean can be regarded as a measure of the 
inner convergence of our photometry. The method ensures 
photometry with no seeing-dependent systematic errors, 
inherent in some other methods, - for images with a PSF 
up to 1'.'4. 

4. Results of VRI photometry 

Our photometry is presented in Tables ^ El and (HI 
The magnitudes were zero-pointed with Yee's (|1988|l ref- 
erence star, with its magnitudes taken from Corrigan et 
al. H1991|l . Our measurements in R band in 1997-2000 are 
plotted in Fig. 2, where the OGLE data (Wozniak et al 
I2000b|l taken in V filter arc shown in grey. For better com- 
parison, our data of Tables 4 and 5 are shifted by small 
amounts along the vertical axis, - 0.1, 0.13, 0.15 and 0.3 
magnitudes for A, B,C and D, respectively. 

The most important brightness changes observed were: 

1. An increase of the image A brightness, starting at the 
end of 1998 and peaking, according to photometry of 
Wozniak et al. lt2000b,l and Alcalde et al. (,20Q2J, in the 
middle of November 1999. We observed almost 0.4™ 
brightening of A image between our observing seasons 
in 1998 and 1999. 

2. A monotonic decline of almost 1.0™ in image B bright- 
ness starting with our earliest, 1995.8 observation. It 
has become the faintest component in R band by 
September 2000. 

3. A strong brightness peak in image C. Our observations 
in July 19-22, 1999 were made near the brightness peak 
of the C component, seen in the well-sampled light 
curves of Wozniak et al. H2000hll . The C image became 
almost 1™ brighter in R band between August 1997 
and July 1999. Thus we have an excellent occasion to 
detect the color change that accompanied the bright- 
ness peak. 

4. A noticeable growth of the D image brightness, which 
is no more the faintest one since September 2000. 

Our measurements, presented in Fig. 2 are in a good 
qualitative agreement with more detailed and accurate 
single-filter light curves of Wozniak et al. I2000bl taken 
in V band for a similar epoch. A large scatter of points 
for 1997 and 1998 in Fig. 2 are due to the dates, when the 
TI 800 X 800 CCD camera was used. We compared our V 
magnitudes, taken with the TI 800 x 800 camera in 1997- 
98, (Table O with the same dates of OGLE monitoring. 



Table 3. A comparison of three programs of Q2237+0305 
photometry: Maidanak (this work, ST-7 CCD), OGLE 
and GLITP; V magnitude differences for A,B,C and D 
components; observations of 2000. 



Programs 


^Va 


AVb 




AVd 


GLITP - OGLE 


0.07 


0.01 


0.14 


-0.18 


Maidanak-OGLE 


0.06 


0.05 


0.07 


-0.15 



and found that the OGLE V magnitudes are systemat- 
ically smaller than our measurements with this camera. 
The greatest differences are for the A and D components, 
reaching approximately 0.2™, with 0.1™ for B and C. As 
seen from comparison of Table 01 with Table O where the 
photometry with the ST-7 camera is presented, the latter 
is almost an order of magnitude more accurate as com- 
pared to the TI 800 x 800 data. 

Our well-sampled and most accurate measurements, 
made in July-September 2000 with the ST-7 camera, - the 
datapoints near the right edge of Fig. ??, - can be seen in 
Fig. 3 more in detail, (see also Table O. For better clar- 
ity, the light curves in Fig. 3 were arbitrarily shifted along 
the magnitude axis, and fit with quadratic polynomials, 
with the 1% error strips shown. Variations of brightness 
of all the components were moderate during this time pe- 
riod, about 0.02™ 0.03™ per month, and may be ap- 
proximated by the second-order polynomials quite well. 
The brightness estimates for the A component are mainly 
within a 1% deviation with respect to the fitted curve. A 
correlation between the rapid brightness variations of all 
the components seen in Fig. 3 could be ascribed to quasar 
intrinsic brightness changes, except that since their am- 
plitudes are larger for the fainter components, these vari- 
ations are probably not real and are more likely due to 
errors. 

We compared our photometry of July-September 2000 
in V band with the OGLE data, obtained for the same 
dates, and, since our data do not overlap with the obser- 
vations of GLITP collaboration, we made a similar com- 
parison between their photometry and that of OGLE. The 
results of such a comparison are presented in Table O 
Here, positive differences mean that OGLE magnitudes 
are smaller. The difference between our photometry and 
that of OGLE program will be even smaller for A, B and 
C images if one takes into account 0.034™ difference in 
magnitude for a star adopted in Wozniak et al. H2000a|l 
and in this work, though systematics for the C component 
will become larger. 

5. Variations of color in Q2237+0305 

The first multicolor observations by Yee l|1988|l have im- 
mediately shown, that the components differ in their col- 
ors. An obvious dependence of the components' reddening 
on the distance to the galaxy nucleus allowed Yee to ex- 
plain it by selective extinction in the dusty matter of the 
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lensing galaxy. This suggestion made it possible to esti- 
mate the extinction law in the lensing galaxy, which, ac- 
cording to Nadeau et al. p991|l and Yee H1988|l . is similar 
to that in our Galaxy. It should be emphasized here, that 
the conclusion was based on the analysis of color differ- 
ences of the components for a fixed epoch. 

As mentioned in the Introduction, a suspicion arose in 
1991 and 1992, that the colors of the components might 
have changed, (Corrigan et al. 119911 Rix et al. 11992(1 . In 
particular, Corrigan et al. did not find any significant vari- 
ations of {B — r) colors of the components with time, 
but they were the first to notice that "there may be a 
small color change in image A as the r magnitude gets 
fainter" (Corrigan et al. 1991^ . They referred to the work 
by Wambsganss & Paczinski 1)1991(1 . where the possibility 
is discussed that, if the quasar structure is wavelength de- 
pendent, microlensing events will differently reveal them- 
selves in different spectral regions. In particular, according 
to Wambsganss & Paczinski fl9915, the bluer inner parts 
of the continuum source might be more strongly amplified 
as compared to the outer parts. 

Rix et al. (19921, analyzing their observations in U 
and R bands with the Hubble Space Telescope, plotted 
their {U — R) colors against {g — i) colors of the compo- 
nents, measured by Yee ((1988(1 . and concluded that they 
"are only marginally consistent" with the reddening line 
derived by Nadeau et al. 1(1991(1 . They suggested, the dis- 
crepancy could be due to either variable dust extinction 
in the lensing galaxy, or to the effects of microlensing 
color changes, first noted by Kayser et al. (1989) and 
later investigated by Wambsganss & Paczynski (jl991.) and 
Wambsganss ((1991)1 in simulations. 

We have already analyzed the behavior of the rela- 
tive colors of the Q2237-I-0305 components qualitatively, 
(Dudinov et al. l20U0al and Dudinov et al. 2000bj), based 
upon our observations on Maidanak in 1995 (Vakulik et 
al. 11997(1 . and in 1997-1998, and also upon all available 
multicolor observations by other authors, (Burud et al. 

0stensen et al. ITMni Rix et al. ITM^ Yee EHHI)- 
A tendency for the components to become more blue as 
their brightness increases has been noted there, but no 
quantitative relationships have been derived. 

We present here our measurements of the colors of the 
A,B,C,D components, and the attempt to quantitatively 
analyze the behavior of {V — R) and {V ~I) color indices of 
the components using our data taken in 1995-2000. VRI 
photometry is presented in Table|Bl and (V — R) and {V — 
I) colors can be seen in Table H Formal errors for these 
quantities, calculated as the errors of the average, range 
from 0.02™ - 0.03" (A component) to 0.03™ - 0.05" (D 
component), for the most accurate observations of 2000, 
and, as seen from Table EJ are within 0.08™ — 0.15™ for 
the observations of 1997-1998, made with the TI 800 x 800 
camera. 

As shown in the previous section, our photometry is 
in quite satisfactory agreement with that obtained by 
other observers for close epochs, (e.g. Alcalde et al. 120021 
Woznyak et al. [2000b|l . At any rate, the discrepancy does 



not exceed that obtained when different algorithms are 
applied to the same data, (Alcalde et al. [2IT02. Burud et 
aL lTOMjl . 

However, one should keep in mind the peculiarities of 
Q2237-)-0305 photometry mentioned above in analyzing 
and interpreting the lightcurves in general, and especially 
those combined from heterogeneous observing data. With 
this in view, more weight should be given to the analysis of 
relative quantities, which are less sensitive to differences 
in observational circumstances and algorithms of image 
processing. In particular, relative colors and relative mag- 
nitudes, as well as their variations are such quantities. 
Examining their behavior in time, and their relationships 
with each other in microlensings can be a valuable source 
of additional information about the physical properties of 
both the quasar and lensing galaxy. In particular, they can 
be used to probe the spatial structure of the quasar at dif- 
ferent wavelengths, (Wambsganss & Paczvnski ll991(l . and 
to determine the extinction law in the lensing galaxy. 

A correlation between {V — I) colors of the compo- 
nents and their R magnitudes can be seen from Fig. 4, 
where the components are marked with different symbols. 
It is interesting to note, that B, C and D components 
are arranged just along a line in this diagram, while the 
A component forms a separate cluster of points. We can 
not refute the possibility of some systematic errors in our 
photometry, but we argue that they would hardly arrange 
the B, C and D components along a single line so well, 
- a correlation reaches 0.8 for them, - and separate the 
A component so significantly. Moreover, we studied the 
systematic errors of our algorithms very carefully in sim- 
ulation and found, that their effect, if present, might only 
slightly bias the color of the C and D components to larger 
values, i.e. make them redder, as compared to A and B, 
(see Vakulik et al. 1997 for more details). We see from Fig. 
4, however, that the cluster of point for A image is shifted 
towards redder colors with respect to the cluster for B, C 
and D. 

If all the components were equally macroamplified, and 
if both the colour differences of the components and their 
variations in time were caused by the interstellar redden- 
ing law, similar to that for our Galaxy, a linear relationship 
between {V — I) and R could be expected, with a regres- 
sion slope of about 0.42 for {V - I) base, (Schild [MZI). 
Microlensing events, with their still unknown brightness- 
color dependence, would disturb and rearrange this order, 
making the components follow the reddening line in the 
average, but forming individual clusters of points, with 
the patterns and stretches, determined by the level of mi- 
crolensing activity at a particular time period, and by the 
unknown character of color-brightness dependence of mi- 
crolensings. 

However, the existing macrolens models predict differ- 
ent macroamplifications for the components. According to 
the macromodel by Schmidt et al. ((1998(1 . rather well con- 
firmed by the observations in emission lines (Fitte & Adam 
[TOMi Racinefl992 Lewis et al. lT^ Saust 1994), and in 
the IR spectral range, (Nadeau, et al. [r991t Agol, et al. 
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121)1)0 1. where no microlensing effects are expected, the A, 
B and D components must be almost equally macroampli- 
fied, with the flux ratios of 0.25, 0.27 and 0.32 respectively. 
The C component is expected to have the least macroam- 
plification factor by this model, - flux ratio of 0.15 is pre- 
dicted for it. It means, that in this case the components 
can not be expected to sit along a line in the {V — I) vrs 
R plot. In the presence of microlensings, the components 
would produce a family of clusters, shifted with respect to 
each other by the amount of flux ratio differences. 

We see in Fig. 4 another situation however. While B, C 
and D components produce three overlapping clusters, all 
of them being stretched approximately in the same direc- 
tion, the A datapoints form a separate cluster, stretched 
along a line with a slope similar to that of the joint B, C 
and D cluster. All the points in the joint cluster are rather 
well correlated, with a correlation index of 0.8 ± 0.1 and 
a regression line slope of 0.33 ± 0.08. The points formed 
by the A component are also rather well correlated, with 
a correlation index of 0.84 and a regression line slope of 
0.36. 

To eliminate possible additive constituents of the color- 
magnitude dependence, which may differ for different com- 
ponents, and to focus on the analysis of changes, we stud- 
ied a correlation between the deviations of {V — I) colors 
from their average over the whole time period, and sim- 
ilarly calculated variations of brightness in R band. The 
diagram can be seen in Fig. 5. The quantities are rather 
well correlated, with a correlation index of 0.75 ± 0.08, 
and a regression slope of 0.31 ± 0.08, in a good agreement 
with that of Fig. 4 for B, C and D components, but the A 
component is not situated separately this time. The un- 
certainties are given for an 80% confidence interval. 

A diagram of color {V — I) - color {V — R), which is 
known to be of great diagnostic importance for the study 
of dust extinction, is usually presented by all authors of 
multicolor observations, e.g., Yec ( 1988), Rix et al. |I992), 
Burud et al. H1998|l . but, as was noted above, only for a 
fixed epoch. In such a diagram, the color indices should 
be proportional to each other for any color base and for 
any type of dust extinction, with the slope determined 
by the reddening law. The diagram, built with the use 
of our measurements (see Table EJ, can be seen in Fig. 
6. A large range of color variations of the C component 
should be particularly noted. It is quite real and can be 
explained by the presence of observations of 1999 in our 
data, - two asterisks near the origin. As was noted, these 
data were obtained near the July 1999 brightness peak of 
C, when it became almost 1™ brighter during two years, 
(Wozniak et al. 12000 bjl . and exceeded the B component in 
brightness. The datapoints in this diagram are found along 
a line with a slope of approximately 1.31 ± 0.14, which is 
much less than 2.15 for these color indices expected for 
the interstellar reddening law in our Galaxy, (Schild|EZI()- 
We conclude that if the extinction law in the lens galaxy is 
similar to that of our Galaxy, the observed color changes 
can not be explained by variable interstellar reddening. 



6. Color Changes Associated with Brightness 
Peaks 

A further perspective of the nature of the observed color 
changes comes from a comparison of the history of bright- 
ness changes with the history of color changes. The general 
features of long-term variations of the components mag- 
nitudes during 1995-2000 in comparison with the simul- 
taneously determined colors can be seen in Fig. 7. Since 
durations of our observing seasons are small as compared 
to the characteristic time scale of the long-term variations 
of colors of the components, we calculated the mean val- 
ues for {V — I) and R for every season, and plotted them 
as a function of time, (the midpoint dates of each season 
are used here). 

Long-term variations of colors of the components are 
clearly seen in this figure, as well as a tendency of the 
components to change their color indices towards smaller 
values (bluer color) as their brightness increases. But the 
tendency is not always straightforward. Some reddening 
of the components, preceding the subsequent decrease of 
their color indices at the stage of component brighten- 
ing can be also seen. As noted above, our observations in 
July 1999 were made very close to the brightness peak of 
C component, that is seen very well from the more com- 
plete and well sampled light curves of the OGLE program, 
Wozniak et al. Ij2000b|) . while the A component was just 
in the middle of the ascending slope of its peak at this 
time, according to the observations of OGLE and GLITP 
programs, (Wozniak et al. I2000bl Alcalde et al lW)^ . 

The relationship between brightness change and color 
is obvious and approximately as expected from models of 
Wanibsganss & Paczynski H1991|> . The most direct corre- 
lation is found for image C, where {V — 1) color is almost 
perfectly anti-correlated with brightness. Thus as the C 
quasar image brightened by almost 1.0™ in R between 
1997.5 and 1999.6, it became bluer by 0.42 magnitudes in 
{V — I) color index (Fig. 7). A second interesting behavior 
is seen in the brightness of image A, where we find that the 
brightness increased by 0.45™ as the color became bluer by 
0.15™ in the {V — I) color index. Just as interesting is the 
color history for image B, which underwent a sustained 
slow brightness drop of 1.0'" during our monitoring pe- 
riod. As it gently declined in brightness, it became 0.25™ 
redder in {V - I) color from 1995.7 to 1998.8, and then 
again became 10% (0.1™) bluer as the brightness contin- 
ued to fade from 1998.8 to 2000.9. 

We have already shown that this is not likely to be pro- 
duced by a hole appearing in some absorbing clouds. If in- 
stead we view the image C brightness peak as a microlens- 
ing artefact where a compact object (star) in the lens 
galaxy passed in front of the quasar and caused the tem- 
porary brightening, we can compare to the calculations 
in Figs 1, 2, and 3 of Wambsganss & Paczinski H1991|l . 
Their models were crafted to apply to Q2237, and they 
show approximately the correct brightness change (1™ in- 
crease in V) and color change (0.4™ bluer in {B — R)) 
for events with 1 year duration, and appear similar to 
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the Wambsganss & Paczinski H1991fl Fig. 2 i,j pattern of 
a quasar image passing outside a cusp of a microlensing 
star. We do not press these calculations further because 
we feel that the failed Wyithe, Turner, and Webster H2000|l 
prediction of a subsequent large brightness change inval- 
idates all models with such simple accretion disc approx- 
imations. However almost any quasar model with an en- 
ergizing central source produces quasar structure which 
is more compact at shorter wavelengths. We expect to 
produce separately a series of models that can reproduce 
the observed effects, based upon the double-ring Schild & 
Vakulik model. 

Although quasar emission lines contaminate the color 
photometry in the continuum-dominated filter bands, we 
doubt that the emission lines are responsible for the large 
brightness-color effects found here, given that the large 
brightness changes observed are always associated with 
microlensing of the quasar continuum. 

7. Conclusions 

1. Our observations demonstrate drastic changes of the 
component magnitudes, which are inherently uncor- 
related in this system, confirming high probability 
for microlensings, predicted for Q2237-I-0305 in 1989, 
(Kayser & Refsdal 11989]) . The highest gradient of 
brightness change was observed for the C component 
between 1997 and 1998, - almost 0.07™ per month in 
our R filter. Almost the same value has been measured 
by OGLE program in their V band for the same time 
period. However, a much more rapid brightness change 
of the C component was detected immediately after 
its extraordinary brightness peak in July 1999 by the 
OGLE program, - almost a 0.2*" decrease per month 
in their V band, (Wozniak et aLEJOOb). 

2. No strong microlensing event occurred in the system 
during our detailed 2.5-months monitoring in July- 
October 2000 but the fact that the B component 
has become the faintest one, after its long continu- 
ous fading beginning in 1995, (see Fig. 2 and TablelSJ. 
No noticeable night-to-night brightness variations were 
detected in this time period. Moderate brightness 
changes were inherent in all the components, reaching 
a 0.03™ decrease per month for B and C, and 0.02™ 
brightening for D, (see Fig. 3 and Table EJ. 

3. All the components demonstrated variations of their 
colors during 1995-2000, which we argue to be real and 
significant. The most prominent change of color was 
observed for the C component, - 0.43™ for its {V — I) 
color index during two years. The {V — I) color indices 
of A, B and D were less variable during the whole 
time period, having changed from 0.3™ to 0.5™ for A, 
from 0.2™ to 0.5™ for B, and from 0.7™ to 0.45™ for 
D component, which became 0.2™ bluer between 1997 
and 2000, having approached the B component in color 
and exceeded it in brightness. 

4. The {V — I), {V — R) color-color plot shown in Fig. ?? 
incorporates all our observations. The regression slope 



is 1.31 ± 0.14 for this diagram, i.e. much smaller than 
a value of 2.16, expected for the reddening line in our 
Galaxy for these color indices. We conclude that the 
brightness and color changes observed are not caused 
by time variations in reddening, but are more probably 
caused by microlensing of source structure that is more 
compact at shorter wavelengths. 

5. The large brightness peak of the C component in July 
1999 was accompanied by large color change in the 
sense that as the C image brightness increased by al- 
most 1™ both in our R and Wozniak et al. (2000b| V 
bands, the color became bluer by 0.43™ in (V—I). This 
is the sense and amplitude expected for microlensing 
of an object that is smaller at shorter wavelengths, and 
modelled previously by Wambsganss and Paczinski 
l|1991|l . The colors of the other components behave 
similarly, though the amplitudes of their color vari- 
ations are smaller, (see Fig. 7). 

6. Returning to Fig. 4, where the relationship between 
the (V—I) colors and R magnitudes of the components 
is shown, we note that the plot is inconsistent with the 
adopted models of macrolensing, e.g. Schmidt et al. 
((T998). We think that most probably the A component 
is macroamplified almost 0.8™ more than B, C and D, 
which have almost equal amplifications. 

We hope, that the data presented here will demon- 
strate the importance of multiband observations of grav- 
itationally lensed quasars in general, and Q2237+0305 in 
particular. More detailed analysis of the obtained data 
and simulation with the new quasar structure model, will 
be presented in the next paper, which is in progress. 
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Figure captions 

Fig. 1. Images of Q2237+0305 for six epochs, obtained in 
R band with the 1.5-m Maidanak telescope; A component 
is at the bottom, C is at the left. 

Fig. 2. Photometry of Q2237+0305 A,B,C,D in R band 

from observations with the 1.5-m Maidanak telescope in 
1997-2000, (large symbols). Photometry in V band by 
OGLE program is also plotted by smaller and fainter sym- 
bols. Our data (Tables 4 and 5) are shifted arbitrarily for 
better comparison, see Sec. 4 for more details. The appar- 
ent brightness discrepancy near the image C 1999 bright- 
ness peak results from the different monitoring filter bands 
used and a significant color change during the brightness 
peak event. 

Fig.3. Photometry of Q2237-I-0305 A,B,C,D in R band 

from observations with the 1.5-m Maidanak telescope in 
2000, July 23 - September 7. The light curves of the com- 
ponents are shifted arbitrarily along the magnitude axis 
with respect to each other for clarity, and approximated 

by second-order polynomials. 

Fig. 4. {V — I) colors vrs R magnitudes, calculated for ob- 
servations of 1995-2000. Note a separate cluster of points 
for A component. A regression line slope for B, C and D 
components is of 0.33 ± 0.08. 

Fig. 5. Variations of (V — I) colors (vertical axis) vrs varia- 
tions of R magnitudes, calculated for observations of 1995- 
2000. A regression line slope is 0.31 ± 0.08. 

Fig. 6. Color {V — I) vrs color {V — R) diagram 
for Q2237-I-0305 A,B,C,D components, plotted from the 
Maidanak data, taken in 1995-2000. A regression line slope 
is 1.31 ±0.14. 

Fig. 7. Long-term variations of R magnitudes (upper 
panel) and {V — I) colors (at the bottom) of A, B, C, 
D components of Q2237-I-0305 from the observations of 
1995-2000. Each point is a result of averaging within one 
observational set. 



Figures 1-7 " are available in "gif format from: 
http: / /arXiv.org/ps/astro-ph/ 
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Table 4. Photometry of Q2237+0305 in R band from observations with the Maidanak 1.5-m telesope in 1995-1999. 



Date 


A 


B 


C 


D 


FWHM(") 


Camera 


Scale("/pix) 


95 09 17 


17 18 ± 03 


17 32 ± 03 


18 13 ± 06 


18 44 ± 07 


0.90 


Pictor 


0.159 


97 07 02 


17 10 ± 05 


17 77 ± 05 

X J . 1 1 — 1 — \j * yj t_/ 


18 10 ± 07 


18 48 ± 10 


0.73 


TI 


0.130 


97 07 03 


1 7 08 + Ofi 


1 7 62 + 06 

X 1 _1_ \J.\J\J 


1 7 98 + 08 


1 8 55 ± 1 2 

A^Kj.tJtJ _1_ Lf.XXj 


0.84 


TI 


0.130 


97 08 29 


17 15 ± 03 


17 72 ± 03 


18 08 ± 12 


18 38 ± 10 


0.85 


TI 


0.130 


97 08 30 


17 14 + 02 


1 7 75 + f) fiQ 


1 o no J_ f) (14 


1041 + 05 


0.71 


TI 


0.130 


97 08 31 


1 7 1 fi + 02 


1 7 7Q + n 03 

X 1 . 1 (_) _1_ KJ.yjtj 


1 Q f)5 + f) nt; 

A.yJ.\JyJ _1_ yj.xjxj 


i Q 40 + n nfi 


0.81 


TI 


0.130 


97 09 01 


17 13 + 04 


17 75 ± 05 

X J ■ 1 — 1 — y^j * yj t^f 


18 01 ± 04 


18.44 ± 0.11 


0.85 


TI 


0.130 


97.10.18 


1 7 20 + 01 


1 7 65 + 02 


1 7 97 + 02 


1 8 39 ± 04 


0.74 


TI 


0.130 


97 10 22 


1 7 2.3 + n 02 


1 7 73 + f) nr. 

X I . 1 fj _1_ \J.\JtJ 


1 Q 04 -L n (14 


1 Q 5(1 + f) 12 


0.84 


TI 


0.130 


97 10 23 


17 18 ± 02 


17 56 ± 04 

X J 1 'J yj — 1 — 


17 96 ± 04 


18 33 ± 06 


0.83 


TI 


0.130 


97.10.24 


17 19 + 03 


1 7 76 + 04 


1 7 96 ± 05 

X I . t l\ J _1_ W . y_l ' J 


1 8 54 ±0 07 


0.76 


TI 


0.130 


97.11.11 


1 7 24 + 02 


1 7 63 + 03 


1 7 95 ± 05 

X I . ti'.i _i_ yj • yi ' J 


18 47 ± 07 


0.77 


TI 


0.268 


97.11.12 


17 25 ± 03 


17 70 ± 05 


17 98 ± 05 

X 1 m\jyj \J *\J%J 


18 46 ± 08 


0.87 


TI 


0.268 


98 07 23 


1 7 1 n + f) 04 


1 7 Q7 + f) f)Q 


1 7 52 + 10 


iQ 22 ± 12 


0.96 


TI 


0.130 


98 07 26 


17 08 ± 04 


17 82 ± 04 


17 43 ± 06 


18 12 ± 07 

X • X ^ — 1 — KJ *\J i 


0.87 


Pictor 


0.159 


98 07 28 


17.14 ± 0.02 


18 00 ± 02 


17 63 ± 06 

X I * yj^_f — 1_ yj w\j yj 


18.32 ± 0.12 


0.88 


TI 


0.130 


98 08 23 


1 7 22 + 07 


1 8 03 + 04 


1 7 41 ± 03 

X ( .^x _i_ yj.yjtj 


18 33 ± 06 


1.00 


TI 


0.130 


98 08 24 


17 18 ± 02 


18 05 ± 04 


17 39 ± 03 


18 30 ± 04 

A^yj yj*\j^ 


1.05 


TI 


0.130 


98 08 25 


17 15 ± 02 


17 97 ± 02 


17 45 ± 02 


18 28 ± 04 

xi^*ijt_7 ^^L- yj*\j^ 


1.08 


TI 


0.130 


98 08 2fi 


17 15 + 003 


1 8 00 + 04 


1 7 45 ± 03 


1 8 38 ± 05 


0.89 


TI 


0.130 


98 08 28 


17 05 ± 02 


17 85 ± 02 


17 45 + 03 

X 1 ♦ i yj _]_ yj 


18 04 ± 04 

A-K^ »yj^ — 1 — yj »yj^ 


0.96 


Pictor 


0.159 


QQ r\Q 9Q 


1 7 n*^ _|_ n no 


1 7 _|_ n no 


1 7 on n no 


1 nn _L n n^ 


u.ou 


r icioi 


n 1 


98.08.30 


17.03 ± 0.02 


17.82 ± 0.02 


17.39 ± 0.03 


17.99 ± 0.04 


0.90 


Pictor 


0.159 


98.08.31 


17.06 ± 0.02 


17.88 ± 0.02 


17.40 ± 0.03 


18.08 ± 0.04 


0.86 


Pictor 


0.159 


98.09.01 


17.08 ±0.02 


17.83 ±0.02 


17.44 ± 0.03 


18.04 ±0.04 


1.04 


Pictor 


0.159 


98.09.02 


17.06 ± 0.02 


17.84 ± 0.02 


17.40 ± 0.03 


18.04 ± 0.04 


0.86 


Pictor 


0.159 


98.10.22 


17.05 ± 0.02 


17.89 ± 0.02 


17.40 ± 0.03 


17.99 ± 0.04 


0.97 


Pictor 


0.159 


98.11.14 


17.07 ± 0.05 


17.89 ± 0.06 


17.43 ± 0.05 


18.07 ± 0.08 


1.05 


TI 


0.268 


99.07.19 


16.76 ± 0.02 


18.01 ± 0.04 


17.11 ±0.02 


18.08 ± 0.03 


0.96 


ST-7 


0.160 


99.07.20 


16.75 ± 0.06 


17.99 ±0.06 


17.14 ± 0.04 


18.15 ±0.11 


1.15 


ST-7 


0.160 


99.07.22 


16.78 ± 0.01 


18.00 ± 0.03 


17.13 ±0.02 


18.09 ± 0.02 


0.87 


ST-7 


0.160 
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Table 5. Photometry of Q2237+0305 in R band from observations in 2000 on Maidanak; 1.5-m telescope, ST-7 CCD 
camera. Seing conditions (FWHM in arcseconds), and number of frames (n) are also presented. 



Date A B C D FWHM n 



00, 


,07, 


.23 


16, 


,729 


± 


0, 


,004 


18, 


,180 


± 


0, 


,009 


17, 


,806 


± 


0, 


,016 


18. 


166 


± 


0, 


.015 


1, 


,08 


16 


00, 


.07, 


.24 


16. 


,726 


± 


0, 


,005 


18, 


.214 


± 


0, 


,011 


17. 


,826 


± 


0, 


,016 


18. 


253 


± 


0, 


.024 


1, 


,36 


28 


00, 


.07, 


.25 


16. 


,731 


± 


0, 


,003 


18, 


.194 


± 


0, 


,008 


17. 


,794 


± 


0, 


,007 


18. 


171 


± 


0, 


.012 


1, 


.11 


30 


00, 


.07, 


.26 


16. 


,734 


± 


0, 


,003 


18, 


.184 


± 


0, 


,009 


17. 


,802 


± 


0, 


,006 


18. 


164 


± 


0, 


.007 


0, 


.78 


8 


00, 


.07, 


.28 


16. 


,743 


± 


0, 


,005 


18, 


.208 


± 


0, 


,005 


17. 


,808 


± 


0, 


,007 


18. 


178 


± 


0, 


.024 


0, 


.89 


7 


00, 


,07, 


.29 


16. 


722 


± 


0, 


,003 


18, 


.175 


± 


0, 


,005 


17. 


799 


± 


0, 


,005 


18. 


163 


± 


0, 


.007 


0, 


.94 


35 


00, 


,07, 


.30 


16. 


,738 


± 


0, 


,003 


18, 


.182 


± 


0, 


,005 


17. 


,790 


± 


0, 


,004 


18. 


184 


± 


0, 


.007 


0, 


.86 


20 


00, 


.08, 


.01 


16. 


,728 


± 


0, 


,010 


18, 


.294 


± 


0, 


,050 


17. 


,887 


± 


0, 


,038 


18. 


328 


± 


0, 


.060 


1, 


.42 


7 


00, 


,08, 


.02 


16. 


,731 


± 


0, 


,003 


18, 


.179 


± 


0, 


,008 


17. 


,788 


± 


0, 


,007 


18. 


152 


± 


0, 


.012 


0, 


.77 


12 


00, 


,08, 


.03 


16. 


,740 


± 


0, 


,004 


18, 


.186 


± 


0, 


,008 


17. 


,819 


± 


0, 


,010 


18. 


171 


± 


0, 


.014 


1, 


.14 


21 


00, 


,08, 


.04 


16. 


,731 


± 


0, 


,004 


18, 


.189 


± 


0, 


,008 


17. 


,803 


± 


0, 


,008 


18. 


171 


± 


0, 


.010 


0, 


.88 


14 


00, 


,08, 


.05 


16. 


,745 


± 


0, 


,002 


18, 


.201 


± 


0, 


,009 


17. 


,832 


± 


0, 


,008 


18. 


188 


± 


0, 


.013 


1, 


.09 


34 


00, 


,08, 


.06 


16. 


,749 


± 


0, 


,005 


18, 


.229 


± 


0, 


Oil 


17. 


,834 


± 


0, 


,007 


18. 


211 


± 


0, 


.013 


1, 


.16 


35 


00, 


,08, 


.07 


16. 


733 


± 


0, 


,008 


18, 


.231 


± 


0, 


,018 


17. 


,867 


± 


0, 


,021 


18. 


266 


± 


0, 


.031 


1, 


.46 


21 


00, 


,08, 


.08 


16. 


726 


± 


0, 


,002 


18, 


.194 


± 


0, 


,010 


17. 


,821 


± 


0, 


,008 


18. 


171 


± 


0, 


.012 


1, 


.07 


23 


00, 


,08, 


.09 


16. 


736 


± 


0, 


,004 


18, 


.188 


± 


0, 


,012 


17. 


,803 


± 


0, 


,009 


18. 


162 


± 


0, 


.010 


0, 


,89 


9 


00, 


,08, 


.10 


16. 


,726 


± 


0, 


,004 


18, 


.182 


± 


0, 


,008 


17. 


,789 


± 


0, 


,010 


18. 


165 


± 


0, 


.010 


0, 


.89 


9 


00, 


,08, 


.11 


16. 


,737 


± 


0, 


,003 


18, 


.193 


± 


0, 


,008 


17. 


,807 


± 


0, 


,012 


18. 


144 


± 


0, 


.009 


0, 


.99 


10 


00, 


,08, 


.12 


16. 


725 


± 


0, 


,004 


18, 


.186 


± 


0, 


,010 


17. 


,800 


± 


0, 


,009 


18. 


143 


± 


0, 


.009 


0, 


,98 


17 


00, 


,08, 


.13 


16. 


,736 


± 


0, 


,002 


18, 


.188 


± 


0, 


,006 


17. 


,811 


± 


0, 


,005 


18. 


139 


± 


0, 


.005 


0, 


,96 


42 


00, 


,08, 


.18 


16. 


,723 


± 


0, 


,006 


18, 


.194 


± 


0, 


,012 


17. 


,807 


± 


0, 


,013 


18. 


136 


± 


0, 


.010 


0, 


,92 


10 


00, 


,08, 


.19 


16. 


,741 


± 


0, 


,004 


18, 


.213 


± 


0, 


,012 


17. 


,857 


± 


0, 


,010 


18. 


154 


± 


0, 


.016 


1, 


.04 


15 


00, 


,08, 


.20 


16. 


,734 


± 


0, 


,003 


18, 


.211 


± 


0, 


,008 


17. 


,809 


± 


0, 


,006 


18. 


167 


± 


0, 


.014 


1, 


.13 


32 


00, 


,08, 


.22 


16. 


,739 


± 


0, 


,004 


18, 


.201 


± 


0, 


,008 


17. 


,825 


± 


0, 


,008 


18. 


158 


± 


0, 


.012 


0, 


,93 


12 


00, 


,08, 


.23 


16. 


,744 


± 


0, 


,005 


18, 


.209 


± 


0, 


,009 


17. 


,820 


± 


0, 


,008 


18. 


165 


± 


0, 


.013 


0, 


.86 


10 


00, 


,08, 


.24 


16. 


,732 


± 


0, 


,003 


18, 


.208 


± 


0, 


,008 


17. 


,828 


± 


0, 


,007 


18. 


167 


± 


0, 


.011 


0, 


.92 


20 


00, 


,08, 


.25 


16. 


,743 


± 


0, 


,005 


18, 


.188 


± 


0, 


,007 


17. 


,825 


± 


0, 


,009 


18. 


135 


± 


0, 


.007 


0, 


.85 


10 


00, 


,08, 


.26 


16. 


,742 


± 


0, 


,008 


18, 


.221 


± 


0, 


,013 


17. 


,831 


± 


0, 


,015 


18. 


176 


± 


0, 


.016 


1, 


.02 


15 


00, 


,08, 


.30 


16. 


747 


± 


0, 


,006 


18, 


.212 


± 


0, 


,010 


17. 


,821 


± 


0, 


,007 


18. 


168 


± 


0, 


.017 


1, 


,12 


17 


00, 


,08, 


.31 


16. 


,736 


± 


0, 


,004 


18, 


.214 


± 


0, 


,007 


17. 


,819 


± 


0, 


,010 


18. 


161 


± 


0, 


.013 


0, 


.96 


12 


00, 


,09, 


.01 


16. 


,737 


± 


0, 


,004 


18, 


.190 


± 


0, 


,009 


17. 


,815 


± 


0, 


,009 


18. 


144 


± 


0, 


.013 


0, 


.85 


20 


00, 


,09, 


.02 


16. 


,752 


± 


0, 


,006 


18, 


.201 


± 


0, 


,012 


17. 


,831 


± 


0, 


,010 


18. 


144 


± 


0, 


.010 


1, 


.00 


10 


00, 


,09, 


.04 


16. 


,749 


± 


0, 


,005 


18, 


.226 


± 


0, 


,017 


17. 


,841 


± 


0, 


,020 


18. 


164 


± 


0, 


.019 


0, 


,92 


9 


00, 


,09, 


.05 


16. 


755 


± 


0, 


,003 


18, 


.220 


± 


0, 


,010 


17. 


,859 


± 


0, 


,010 


18. 


152 


± 


0, 


.013 


1, 


,16 


21 


00, 


,09, 


.06 


16. 


,757 


± 


0, 


,003 


18, 


.245 


± 


0, 


,011 


17. 


,863 


± 


0, 


,012 


18. 


162 


± 


0, 


.012 


1, 


,13 


12 


00, 


,09, 


.07 


16. 


,746 


± 


0, 


,005 


18, 


.207 


± 


0, 


,006 


17. 


,833 


± 


0, 


,006 


18. 


120 


± 


0, 


.008 


0, 


,80 


10 


00, 


,09, 


.08 


16. 


,743 


± 


0, 


,004 


18, 


.204 


± 


0, 


,010 


17. 


,837 


± 


0, 


,011 


18. 


140 


± 


0, 


.012 


1, 


,13 


22 


00, 


,09, 


.09 


16. 


,744 


± 


0, 


,004 


18, 


.223 


± 


0, 


,012 


17. 


,834 


± 


0, 


,008 


18. 


142 


± 


0, 


.021 


0, 


,97 


10 


00, 


,09, 


.10 


16. 


,762 


± 


0, 


,003 


18, 


.235 


± 


0, 


,026 


17. 


,864 


± 


0, 


,013 


18. 


175 


± 


0, 


.012 


0, 


,96 


10 


00, 


,09, 


.23 


16. 


,762 


± 


0, 


,005 


18, 


.238 


± 


0, 


,018 


17. 


,849 


± 


0, 


,012 


18. 


095 


± 


0, 


.027 


1, 


,04 


11 


00, 


,09, 


.24 


16. 


,746 


± 


0, 


,010 


18, 


.264 


± 


0, 


,028 


17. 


,851 


± 


0, 


,026 


18. 


151 


± 


0, 


.027 


1, 


,06 


10 


00, 


,09, 


.26 


16. 


,772 


± 


0, 


,005 


18, 


.253 


± 


0, 


,017 


17. 


,855 


± 


0, 


,009 


18. 


118 


± 


0, 


.018 


1, 


,16 


10 


00, 


,09, 


.27 


16. 


,749 


± 


0, 


,005 


18, 


.240 


± 


0, 


,008 


17. 


,845 


± 


0, 


,008 


18. 


111 


± 


0, 


.013 


1, 


,18 


20 


00, 


,09, 


.28 


16. 


,764 


± 


0, 


,005 


18, 


.273 


± 


0, 


,015 


17. 


,845 


± 


0, 


,010 


18. 


143 


± 


0, 


.017 


1, 


,04 


10 


00, 


,09, 


.30 


16. 


,784 


± 


0, 


,008 


18, 


.265 


± 


0, 


,019 


17. 


,909 


± 


0, 


,020 


18. 


149 


± 


0, 


.021 


1, 


,21 


10 


00, 


,10, 


.07 


16. 


,776 


± 


0, 


,012 


18, 


.294 


± 


0, 


,037 


18. 


,001 


± 


0, 


,034 


18. 


075 


± 


0, 


.033 


1, 


,52 
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Table 6. VRI Photometry of Q2237+0305 in 1995-2000; Maidanak, 1.5-m telescope. 



Date 


A 


B 


C 


D 


Sp.band 




17.34 ± 0.04 


17.44 ± 0.03 


18.41 ±0.10 


18.66 ± 0.08 


V 


95.09.17 


17.18 ±0.03 


17.32 ±0.03 


18.13 ±0.06 


18.44 ±0.07 


R 




16.99 ±0.03 


17.21 ±0.04 


17.83 ± 0.08 


18.23 ±0.09 


I 




17.43 ± 0.03 


17.95 ±0.05 


18.42 ±0.05 


18.78 ±0.07 


V 


97.08.29 


17.15 ±0.03 


17.72 ± 0.03 


18.08 ±0.12 


18.38 ±0.10 


R 




16.96 ± 0.03 


17.60 ± 0.04 


17.81 ± 0.04 


18.15 ±0.07 


I 




17.40 ± 0.03 


17.94 ±0.04 


18.35 ±0.05 


18.79 ±0.06 


V 


97.08.30 


17.14 ±0.02 


17.75 ± 0.03 


18.08 ± 0.04 


18.41 ± 0.05 


R 




16.89 ± 0.02 


17.55 ± 0.04 


17.82 ± 0.03 


18.11 ±0.06 


I 




17.47 ±0.05 


18.01 ± 0.07 


18.38 ± 0.06 


18.91 ±0.09 


V 


97.09.01 


17.13 ±0.04 


17.75 ± 0.05 


18.01 ± 0.04 


18.44 ±0.11 


R 




16.93 ±0.03 


17.61 ± 0.04 


17.83 ±0.04 


18.19 ±0.09 


I 




17.30 ±0.02 


18.11 ±0.05 


17.66 ± 0.05 


18.52 ± 0.05 


V 


98.07.26 


17.08 ± 0.04 


17.82 ± 0.04 


17.43 ± 0.06 


18.12 ±0.07 


R 




16.92 ±0.03 


17.71 ±0.04 


17.32 ±0.04 


17.87 ±0.06 


I 




17.44 ± 0.03 


18.23 ±0.03 


17.78 ±0.03 


18.64 ±0.05 


V 


98.07.28 


17.14 ±0.02 


18.00 ±0.02 


17.63 ±0.06 


18.32 ±0.12 


R 




16.93 ± 0.04 


17.72 ±0.05 


17.47 ± 0.05 


17.82 ± 0.06 


I 




17.52 ±0.02 


18.23 ± 0.02 


17.77 ± 0.04 


18.84 ± 0.07 


V 


98.08.23 


17.22 ±0.07 


18.03 ±0.04 


17.41 ±0.03 


18.33 ±0.06 


R 




17.09 ± 0.04 


17.83 ± 0.03 


17.36 ± 0.04 


18.14 ±0.20 


I 




17.34 ± 0.05 


18.28 ± 0.05 


17.65 ± 0.06 


18.48 ±0.11 


V 


98.11.14 


17.07 ± 0.05 


17.89 ± 0.06 


17.43 ± 0.05 


18.07 ±0.08 


R 




16.94 ± 0.03 


17.79 ± 0.05 


17.44 ± 0.03 


17.92 ± 0.06 


I 




16.89 ±0.04 


18.12 ±0.05 


17.19 ±0.03 


18.34 ±0.11 


V 


99.07.20 


16.75 ±0.06 


17.99 ± 0.06 


17.14 ±0.04 


18.15 ±0.11 


R 




16.62 ± 0.02 


17.76 ± 0.05 


17.07 ± 0.02 


17.84 ± 0.04 


I 




16.92 ±0.02 


18.16 ±0.03 


17.22 ±0.02 


18.40 ±0.08 


V 


99.07.22 


16.78 ±0.01 


18.00 ±0.03 


17.13 ±0.02 


18.09 ±0.02 


R 




16.63 ±0.01 


17.79 ± 0.03 


17.09 ± 0.02 


17.82 ±0.03 


I 




16.888 ± 0.005 


18.401 ± 0.016 


18.011 ± 0.008 


18.400 ± 0.007 


V 


00.07.26 


16.734 ± 0.003 


18.184 ± 0.009 


17.802 ± 0.006 


18.164 ± 0.007 


R 




16.577 ± 0.004 


17.973 ± 0.007 


17.590 ±0.008 


17.911 ±0.006 


I 




16.880 ± 0.004 


18.387 ±0.015 


17.992 ±0.011 


18.392 ±0.017 


V 


00.08.04 


16.731 ± 0.004 


18.189 ± 0.008 


17.803 ± 0.008 


18.171 ±0.010 


R 




16.577 ±0.004 


17.968 ±0.007 


17.604 ± 0.006 


17.919 ±0.004 


I 




16.894 ± 0.008 


18.423 ±0.010 


18.032 ±0.017 


18.381 ±0.020 


V 


00.08.09 


16.736 ± 0.004 


18.188 ±0.012 


17.803 ± 0.009 


18.162 ±0.010 


R 




16.586 ± 0.005 


17.976 ± 0.008 


17.616 ±0.010 


17.925 ±0.016 


I 




16.878 ± 0.008 


18.365 ± 0.025 


18.021 ± 0.030 


18.350 ± 0.025 


V 


00.08.18 


16.723 ±0.006 


18.194 ±0.012 


17.807 ±0.013 


18.136 ±0.010 


R 




16.574 ±0.006 


17.975 ±0.010 


17.612 ±0.013 


17.899 ±0.012 


I 




16.882 ± 0.006 


18.369 ± 0.014 


18.008 ± 0.010 


18.370 ± 0.015 


V 


00.08.25 


16.743 ± 0.005 


18.188 ± 0.007 


17.825 ± 0.009 


18.135 ±0.007 


R 




16.578 ± 0.003 


17.982 ± 0.008 


17.616 ± 0.004 


17.884 ±0.011 


I 




16.886 ± 0.006 


18.398 ±0.013 


17.997 ±0.007 


18.328 ±0.016 


V 


00.09.07 


16.746 ± 0.005 


18.207 ± 0.006 


17.833 ± 0.006 


18.120 ±0.008 


R 




16.583 ± 0.006 


17.988 ± 0.010 


17.612 ± 0.008 


17.884 ± 0.012 


I 




16.918 ±0.005 


18.443 ±0.011 


18.057 ±0.013 


18..365 ± 0.016 


V 


00.09.28 


16.764 ± 0.005 


18.273 ± 0.015 


17.845 ± 0.010 


18.143 ±0.017 


R 




!().()()() ± ().()()() 


18. ()()() ± 


17.()57 ± ().()()() 


17.871 ± O.OK) 


1 
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Table 7. {V - R) and {R - I) colors of Q2237+0305 A,B,C,D; 1995-2000, Maidanak, 1.5-m telescope. 



Date 


A 




B 






C 


D 
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V-I 
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u.oo 


A7 
U.^ / 


U. / z 


9fi 07 QS 
zo.u ( .yo 


99 
u.zz 


'^R 
U.oo 


9Q 

u.zy 


/10 
U.4U 


9*^ 
U.ZO 


'^A 
U.O^ 


AC\ 
U.^U 


u.oo 


9S 07 QS 

Zo.u / .yo 


"^0 


p;i 

U.Ol 


9^ 
U.Zo 


p;i 

U.Ol 


1 
u.io 


"^l 
U.Ol 


"^9 
U.OZ 


n 70 

U. ( u 


9^ OS QS 


O QO 


O /1Q 

U.'iO 


O 90 

u.zu 


n ^^o 
u.ou 


u.oo 


O A'l 


O Cil 
U.Ol 


n 7n 
u. r u 


1/1 11 QS 
1^.1 1 .yo 


97 
U.Z / 


/lO 


u.oy 


^0 

u.oy 


99 
u.zz 


91 
U.Z 1 


0/11 
U.^1 


u.oo 


20.07.99 


0.14 


0.27 


0.13 


0.34 


0.05 


0.12 


0.19 


0.50 


22.07.99 


0.14 


0.29 


0.16 


0.37 


0.09 


0.13 


0.31 


0.58 


26.07.00 


0.15 


0.31 


0.22 


0.43 


0.21 


0.42 


0.24 


0.49 


04.08.00 


0.15 


0.30 


0.20 


0.42 


0.19 


0.39 


0.22 


0.47 


09.08.00 


0.15 


0.31 


0.23 


0.45 


0.23 


0.42 


0.22 


0.46 


18.08.00 


0.15 


0.30 


0.17 


0.39 


0.21 


0.41 


0.21 


0.45 


25.08.00 


0.14 


0.30 


0.18 


0.39 


0.18 


0.39 


0.23 


0.49 


07.09.00 


0.14 


0.30 


0.19 


0.41 


0.16 


0.38 


0.22 


0.44 


28.09.00 


0.15 


0.32 


0.17 


0.44 


0.21 


0.40 


0.22 


0.39 
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